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Abstract 

We investigate the discrete chiral transformation of a Majorana fermion 
on a torus. Depending on the boundary conditions the integration measure 
can change sign. Taking this anomalous behavior into account we define a 
chiral order parameter as a ratio of partition functions with differing bound- 
ary conditions. Then the lattice realization of the Gross-Neveu model with 
Wilson fermions is simulated using the recent 'worm' technique on the loop 
gas or all-order hopping representation of the fermions. An algorithm is 
formulated that includes the Gross-Neveu interaction for N fermion species. 
The critical line rndg) is constructed for a range of couplings at = 6 and 
for N = 2, the Thirring model, as examples. 
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1 Introduction 



Many textbooks on quantum field theory start with scalar fields and discuss, after 
preparing the arena with non-interacting free fields, the quartic self-coupling as 
the simplest conceivable interaction. At this stage renormalization or, on the 
lattice, the continuum limit has to be discussed in a non-trivial way. A mild 
generalization leads to real scalar fields self-coupled by a unique quartic O(A^) 
invariant term. As is well known the kinetic term dictates the field dimension and 
implies a mass dimension 4 — D for the coupling constant. Thus the self-coupled 
scalar theory is renormalizable in D = 4. In fact, D = N = Ais the standard model 
in the limit where all fields except the Higgs are stripped away. The perturbative 
renormalization group /5-function is positive for small coupling and the models 
are hence perturbatively trivial. Simulations with simple lattice discretizations 
suggest that this features is also true nonperturbatively. 

An analogous fermionic model can be written in terms of A^ 'real' Majorana 
fermions. Now, because there is only one derivative in the fermionic kinetic term, 
a quartic coupling has dimension 2 — D. The renormalizable theory in D = 2 
has in fact again a unique O(A^) invariant interaction. It has been proposed and 
considered in the limit of A^ cxo a long time ago [Ij and is now known as 
the Gross-Neveu (GN) model. For the sign of the quartic interaction that has 
been discussed to be the physically meaningful one in [T] the GN model is even 
asymptotically free. This feature that is shared by QCD was one of the main 
reasons for Gross and Neveu to study the model as an analog toy-model. The 
analogy goes even further as there also is a chiral symmetry in the GN model. 
Although it is only discrete, this symmetry involves the two-dimensional analog 
of 75 and is only present for massless fermions, which are thus distinguished by 
an enhanced symmetry. Gross and Neveu have shown that the chiral symmetry 
breaks spontaneously (at infinite A^) and that mass generation and 'dimensional 
transmutation' take place. This beautiful parameter-free theory is then an analog 
of QCD in the chiral limit. 

There have been some early efforts [2j to study this model beyond perturba- 
tion theory by Monte Carlo simulation. Attempts employing today's standard 
techniques used for QCD ^ were made in |1] and [Sj. However in [6J and [7j we 
have proposed techniques to simulate the fermionic Grassmann integral in a more 
direct and efficient way using its loop gas representation [H], see also P] for a related 
effort. In this formulation, using Wilson fermions, the simulations resemble those 
for simple spin models with efficient algorithms. No nearly singular inversions and 
no nonlocal effective bosonic actions occur. While in more than two dimensions a 
sign problem still prevents the application of this technique [7] it is well suited for 
the GN model as we demonstrate below. 

Like in QCD, the important chiral symmetry is violated by the Wilson dis- 
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cretization and an additive mass renormalization, i.e. a nontrivial definition of 
tlie critical point, is required. In tliis article we discuss, that, depending on the 
boundary conditions on a finite torus, an extra minus sign from the non-invariance 
of the Majorana fermion measure has to be taken into account in symmetry rela- 
tions. We consider this as (a primitive form of) an anomaly. One such symmetry 
relation allows us to define a finite volume chiral order parameter. This way of 
defining the chiral point in parameter space is particularly convenient in the loop 
gas representation, where it corresponds to a well measurable topological variable. 

After a definition of the model in section 2 the anomaly is discussed for the 
continuum theory in section 3 and connected with the fluctuating boundary con- 
ditions in section 4. In section 5 the loop gas representation and simulation on 
the lattice are described and numerical results are reported in section 6. Finally 
we end in section 7 with conclusions and an outlook. 

2 O(A^) invariant Gross Neveu model 

The O(A^) invariant GN model [1] is best written in terms of self-coupled Ma- 
jorana fermions C,a,a{x) with the Euclidean action 

S = J d^ai,d, + m)^-^jm^ (1) 

where the contracted index a = 1, . . . ,N refers to the global 'flavor' symmetry 
and a = 1, 2 is spin. For a Majorana fermion ^ is not independent but just an 
abbreviation for 

^=CC (2) 

with the 2-dimensional Dirac matrices 70, 7i and the charge conjugation matrix C 
satisfying 

jJ = -C^^C-\ C = -C^, {7^,7.} = 2V, il,y = l,- (3) 

In two dimensions the model is renormalizable in the strict sense, no other four- 
fermion interaction is compatible with the global O(A^). Moreover, the models 
with N ^ 3 are perturbatively asymptotically free. The Thirring model N = 2 
has no divergent coupling renormalization. At vanishing fermion mass m = 0, 
there is in addition a discrete chiral invariance 

^ l5^, 75 = no7i, ? ^ -hb]- (4) 

A standard step in treating the GN model, both in the continuum and on the 
lattice, is to factorize the interaction into bilinears at the expense of introducing a 
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scalar a{x). Foreclosing for a moment the lattice version, where d'^x becomes the 
cell volume a^, we write in the Boltzmann factor at each x 



^^mr = r ^^-4-'-49<i_ (5) 
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This may be formally proven by shifting a, but perhaps more convincingly, we 
state that both sides agree as polynomials in the nilpotent Grassmann bilinear 

After the introduction of cr{x) we are led to the fermionic partition function 

Z{m + ga; e) = j D.^e'-^ ^ 'i'^~ii^^^d,+m+ga)ti ^ pf [^(7^9^ + m + ga)] (6) 

as the basic fermionic building block. In this formula ^ represents only one flavor. 
We imagine this fermionic functional integral here formally on a finite torus of 
extent Lq x Li in the continuum. The two component ('bits') object e with G 
{0, 1} labels four possible (anti)periodic boundary conditions 

ax±fiL^) = {-lY-a^). (7) 

The Pfaffian Pf is the counterpart of the determinant associated with Dirac fermions. 
Note that Pf is only defined for antisymmetric matrices and that the operator in 
([6]) [absorbing m and g now into cr] 

A(a)=C(7A + ^) = -^^M (8) 

indeed has this property for all four choices of e. One relation for Pfaffians is that 
they square to determinants, 

(Pf[A(a)])^ = Det[A{a)]. (9) 

Note that the sign of Pf , which is our focus in the next section, is not determined 
by the determinant. 



3 Symmetry of the Pfaffian in a finite volume 

By a suitable choice of Dirac matrices A can be taken to be real. Hence its 
eigenvalues are imaginary and come in pairs 

A{a)f = tXf, A{a)r = -i\r. (10) 

In addition, by using 75, it is trivial to demonstrate that A{—a) possesses precisely 
the same pairs of eigenvalues. By arguments analogous to those used for Weyl 
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fermion^ in connection with the global anomaly [10] a functional definition of the 
Pfaffian as a root of the determinant has to be smooth under deformations of a. 
For some reference field a where there are no vanishing eigenvalues we define the 
Pfaffian by multiplying all positive eigenvalues A only. This represents a certain 
definition of the (irrelevant) over-all phase. A possible simple reference field would 
be (j{x) = m > 0, the free massive Majorana fermion. Then the Pfaffian for other 
field configurations including its sign follows from smoothness. 

Again as in this requires us to follow the fiow of individual eigenvalues as 
functionals of a. As a pair of them crosses zero the associated eigenvectors allow 
us to follow the members of the pair separately and defines a unique continuation 
of the phase of the Pfaffian through the crossing. From the spectral properties dis- 
cussed above we conclude that Z[a\ e\ must be either symmetric or antisymmetric 
in a 

Z{a- e) = Ti{e)Z{-a- e), 7r(£) e {+1, -1}. (11) 
In the next paragraph we shall show that 

1 +1 for £^ = (1,0), (0,1), (1,1) 



holds. 

We consider an interpolation 

g{t) = Z{ta; e) = Ff[A{ta)], t G [-1, 1] (13) 

and assume Z[a;e] 7^ as otherwise the proof is trivial. Then the question of 
(anti) symmetry in t can be 'decided' in an infinitesimal neighborhood of t = 
which we can control by leading order perturbation theory. 

Except for = (0, 0) the eigenvalues at t = do not vanish but are of order 
L~^. Hence there is no zero-crossing, which proves the lower line of (|T2l) . In the 
totally periodic case there is a pair of x-independent zero-modes. The 2-spinors 
u± are chosen to diagonalize the spin matrix of the perturbation Cta 

Cu± = ±m±, A{0)u± = 0. (14) 

Now in leading order degenerate perturbation theory we find two eigenvalues 

t /• . 



A± ~ ±— — / d'xa{x) (15) 

associated with these eigenvectors. 

^We admit here that our arguments about infinite objects disregarding renormalization are 
at a highly formal level. However, this will be substantiated on the lattice in our case. 
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Unless this integral vanishes we see that exactly one pair of eigenvalues crosses 
zero as we move t = + e— >-t = — e and hence the Pfaffian changes sign. The 
sign TT cannot jump with a. Thus, if the integral above vanishes, i.e. a has a 
vanishing zero- momentum component, we deform a a little such that this is not 
the case anymore. Now fll2l) is established. 

We note that our general result above is consistent with the explicit evaluation 
in the continuum limit of the Wilson-Pfaffian for constant a, see appendix A in [6J. 
Inspired by this calculation we give another formal proof of the above result for a 
continuum theory regularized by a cut-off in momentum space. We start from the 
identities 

75^A(a)75 = A[-a) (16) 

and 

Pf[75^A(a)75] = Det75Pf[^(fx)]. (17) 

The Det in the last line is taken in the function space over which ^q,(x) is integrated 
(with anticommuting expansion coefficients). If Pf is written as a Grassmann 
integral, then Det 75 is the Jacobian of ^ — >■ 75^. Each 2x2 block contributes a 
factor det 75 = — 1 in the usual sense. We now imagine to regularize the fermion 
integral by a finite mode number cutoff. If we impose this cutoff in a way preserving 
the symmetry of rotations by vr then for each incorporated mode also — has to 
be included. This leads to an odd total mode number only for = (0, 0) because 
of the unpaired mode = (0, 0). Thus we have 

Det75 = 7r(£) (18) 

representing a non-invariance of the integration measure in the case of fully periodic 
boundary conditions. It also follows from the foregoing discussion that expectation 
values in the GN model at odd with such boundary conditions are ill defined at 
the chiral point in any finite volume, because the normalizing partition function 
vanishes. This arises at least superficially in a similar fashion as with the global 
anomaly in [TU] . 

4 Observables from fluctuating boundary condi- 
tions 

4.1 Chiral order parameter 

We have introduced four different boundary conditions for the fermions in the 
GN model. We may consider the choice among them as an additional dynamical 
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variable and consider general superpositions with an amplitude z{e) 

Z{a;z) = Y,4e)Zia;e). (19) 
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Below we shall be interested in two special choices for the amphtudes 

= ^ - 4,(0,0), z^{e) = ^. (20) 

The first choice is distinguished by leading to a strictly positive hopping or world- 
line expansion for the finite volume GN model on the lattice [B]. The anomalous 
symmetry (including the sign of Det75) found in the previous section implies the 
relation 

Zi-a;z+)=Zia;z^). (21) 

For either choice — jointly for all flavors — we find partition functions for the 
interacting GN model 

Z2''{m,g) = J Dae-y'^"'''"[Z{m + ga;z±)f (22) 

which coincide at the chiral point m = as one shows by transforming a ^ —a. 
We hence may define an order parameter that vanishes at the chiral point 

X = j^HZ^Vzr)- (23) 

More conventional order parameters may be based on chirally odd observables 
^{x)^{y). In view of our later simulations we define 

{Uiu)-^,,iv)). = ^ j Dae-y^^\z{m + ga-z)f-' (24) 
X ^z{e) I D,^e-^^'''"«(>^''+"+5'^)«ea(M)e/3(^^)- 

Here we indicated the periodicity in the integration measure. Note that we always 
normalize with respect to the z+ ensemble. Such correlations are 2L^-periodic 
but neither periodic nor antiperiodic over the original torus while bilinear local 
densities remain periodic. Chiral symmetry implies 

75(^aa(^i)|^(t^))^±75 = - {^aa{u)l p^,{v)) (25) 

An integration over a 2Li long interval projects to the spatially periodic component 
at zero spatial momentum and we obtain the chirally odd scalar ks 
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with both ^ indices contracted here. For the purpose of monitoring the field 
normahzation we define in addition 



1 



+Li 



kvixo) = -- dx,{aOhoax))z+- (27) 

We only use z+ boundary conditions here, as they will lead to the highest precision 
on the lattice. In this way ky has both chirally even and odd components. Due to 
parity, ks is an even function of xq while ky is odd. 

We form the ratio ks/ky to cancel multiplicative renormalization factors. While 
with exact chiral symmetry ks vanishes, with a regularization that violates chi- 
ral symmetry by effects linear in the lattice spacing a we can only expect to find 
ksi^o) = 0(a) if X = has been imposed as a renormalization condition. We 
therefore tie Xq to a physical length scale to obtain a scaling situation and define 

ks{nLo/4:) 

MV^ (- = 1,2,3) (28) 

as universal quantities vanishing in the continuum limit at the chiral point. The 
finite size Lq supplies an infrared regulator and keeps our model well-defined (for 
the boundary conditions chosen) at the chiral point in a finite volume. This is sim- 
ilar to the Schrodinger functional of QCD [llj, where the PCAC relation replaces 
the discrete chiral symmetry relations invoked here. 



5 Lattice, Wilson fermions, loop gas 

To pass to the lattice formulation we restrict ^ (x) to the sites of a toroidal square 
lattice with spacing a and integer (anti)periodicity lengths L^/a. In addition we 
replace the Dirac operator in and ([H]) by the Wilson matrix (with r = 1) 

l,d, 7m^m - ^d,d;. (29) 

We use here standard notation and d^,d*^,d^ denote the forward, backward and 
symmetrized nearest neighbor difference operators. The second term eliminates the 
doubler modes but also breaks the chiral symmetry (jll) just like a mass term. The 
consequences of chiral symmetry however are expected to emerge in the continuum 
limit. To reach it one has to determine a critical value mc{g) as there is no 
symmetry to prevent an additive mass renormalization. From here on we use 
lattice units and set a = 1. 
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5.1 Loop gas representation of the GN model 

The Wilson-discretized GN model has aheady been studied in [5] by hybrid 
Monte Carlo methods and Major ana- Wilson fermions were simulated in [6] by a 
cluster algorithm in the loop gas representation elaborating on [8]. Most recently 
we have studied a worm algorithm ^7] and we will now make extended use of results 
from this paper which should be consulted for details. 

In [6] we have shown that the Grassmann integral (IMl) — for one flavor at this 
point — can be written as a sum over contributions from non-intersecting loops 
on the lattice with one additional open string connecting the locations of the field 
insertions u and v. The underlying loop ensemble is 

Z= Q{k;u,v)2-^/^ Yl {2 + m + g(T{x)). (30) 

u,v,{k{l)} xeM[k] 

In this formula the link variables k{l) = 0, 1 are one on the loops/string and zero 
elsewhere. We also say that links with k{l) = 1 carry dimers. The constraint 
Q{k; u,v) = 0, 1 restricts the k configurations to the subset corresponding to the 
hopping expansion of the Majorana Wilson fermion with field insertions at u,v. 
It demands that zero or two dimers are attached to all sites except u, v if they do 
not coincide. In this case they must see exactly one dimer each, the beginning and 
end of the string connecting them. The integer C counts corners, sites connected 
to two orthogonal dimers. Finally the monomer set Ai [k] for each configuration k 
consists of all sites with no dimers attached. 

For each /c-configuration we define topological variables e^[/c] G {0, 1} as fol- 
lows. For each direction fi on the torus we distinguish one sheet orthogonal to the 
yU-direction (a line for D = 2) of links in that direction. This is the same construc- 
tion that is used in connection with the center symmetry in lattice gauge theory. 
A definite convention is given in [7j. Then is zero or one depending on whether 
an even or odd number of dimers belong to the sheet. Thus we decompose the loop 
gas into four classes, see pi- If a contribution has u = v,9{k,u,u) = 1 there are 
only closed loops and coincides with their total mod 2 winding number in the 
respective direction. We also associate possible antiperiodic boundary conditions 
with minus signs on bonds in the respective sheets where we 'close the torus'. Now 
the fermionic partition function is given by 

Z{a;z) = J2<dik- 0,0)2-^/' J] (2 + cT(a;))$(e; z) (31) 

k x€Mlk] 

where we took u = v = to only sum over closed loop 'vacuum' configurations. 
The amplitude $ is 

$(e; z) = (24,(0,0) - (32) 
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The first factor stems from spin and Fermi statistics in the Grassmann integral 
(see [6], [7]). In particular it is +1 for all topologically trivial loop configurations 
which amounts to the absence of a 'hard' sign problem in two dimensions. The 
second factor corresponds to the fluctuating boundary conditions. If we now insert 
z = z± from (1201) we find 

$(e;z+) = l, <l>(e;2;_) = 25,,(o,o), (33) 

and in particular the announced positivity for the 2+ ensemble. 

To pass to the partition function of the interacting GN model we insert the 
loop representation simultaneously for all flavors and find 

N 

Zf{m,g)= J2 l[^Q{ka,0,0)2~^^/'<l>{ea,z^)]l[c{M{x)). (34) 

In this formula ka, Ca, are defined as before, just for each flavor independently. 
The field M{x) with possible values 0, 1, . . . , depends on all ka (left implicit) 
and yields the total number of monomers at x, i.e. it counts to how many A^[fca] 
the site x belongs. The various flavors thus interact through the factor 

-^e-'^'''[2 + m + ga]^. (35) 
-oo v27r 

The weights c(M) may conveniently be generated from a recursion formula. In 
the Monte Carlo updates below we shall rather need ratios 

^(-)^^ <-) 

which are directly generated by 

r(0) = 2 + m, r(M) = 2 + m + — f- — -. (37) 

r[M — \) 

5.2 GN worm algorithm 

Finally in our simulation we sample the ensemble 

^GN^ ^ e(A;,;w,t;)2-^^/2j][e(fc,;0,0)2-^"/^]nc(M(x)). (38) 

The summation now includes a sum over j = 1, 2, . . . , which is an 'active' flavor 
to which the insertions u, v refer. A valid update procedure is given by simply 
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using the algorithm of [7j for the active flavor with the external field replaced 
c(M(x)) which depends also on the non-active flavors that are momentarily frozen 
during the update. Whenever u = v is encountered the active flavor is changed 
randomly and ergodicity is achieved. We summarize the whole procedure in the 
following paragraph. 

The simulation alternates between analogous w-updates and f -updates, of which 
we only detail the latter. In addition we alternate between type I and II. Type I 
involves: 

• Pick a nearest neighbor v' of v, call the link / = {vv') 

• Propose the combined move v v', kj(l) — 1 — kj{r) 

• Compute the ratio q of the weight given in [7j after and before the move where 
we suitably substitute r(M) for the external field ratio. The Metropolis 
decision is controlled by 

In step II one only acts if m = f is met. In this case we: 

• Pick a random new value j G {1,2,..., A^} for the active flavor. We may 
actually imagine this as a Metropolis proposal tried after every type I move. 
Whenever u ^ v holds it will be rejected due to the ^-constraints while for 
M = f it is accepted with probability one. 

• With probability pkick "we re-locate u = v randomly on the lattice. 



5.3 Observables in the loop gas representation 

We define expectation values in the ensemble (|38ll 

m) = ^ E OQ{kf,uM'^-^^''\\[Q{K-.m2-^^'']^^^ (39) 

i,u,v,{ka} ai^i a: 

where O may depend on all quantities that are summed over. The chiral order 
parameter x pSl) is easily measured by 

In other words, at the chiral point x = the fraction of all vacuum [u = v) 
configurations has the trivial topology in all flavours simultaneously. In the trivial 
case (7 = this arises for m = from independent factors 1/2 for each species. 
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Other observables considered in this paper are the correlations 



{ia{x)U{^))z = Sab TTT-^T • (41) 

Here 5^^'^ is the £-(anti)periodic lattice 5 function and the spin matrij^iS has been 
defined in [7]. It carries two spin indices that we suppressed at The sign 

tpj collects the phases from spin matrix elements and fermion loops in the active 
flavor, see eq. (48) in [7j. It is evolved during the updates and thus available at 
every step. Applying this correspondence for ks and ky in ( l26l) . ( 1271) . using ( 133|) . 
we arrive at the translations 

, . . ((^,tr[7o5(A;,; ..,.;)] [(-l)-o>^4o!.o-.o - 41.-.o])) .... 



and 



with 



, , , ((V^, tr[5(fc,;^,t;)][(-l)-°>^4o!.o-.oi^; - Cq-.q^^/])) ,,o, 
ks{xo) = (43) 



i?f = -(l±[|[24.,(o,o)]l. (44) 

At Xq = Lq/2, due to parity covariance, S^^^ in ky and S^^^ in ks do not contribute. 
Although these expressions may look somewhat intimidating, all quantities enter- 
ing are readily available from the updates and the accumulation proceeds quickly 
in the form of histograms with bins for each time separation uq — vq and each of 
the four possible Cj. To debug the code we have found it useful to run at 5^ = 
comparing with the known exact answers. 



6 Numerical tests in the Gross-Neveu model 

We have generalized the code described in [7] to include an arbitrary number of 
flavors coupled by the interaction term flH3|) and we thus sample the ensemble (EH])- 
We have accumulated the necessary observables entering into x and Xn,n = 1,2,3. 
By consulting perturbation theory and some preliminary experiments, we have run 
series of simulations varying L at fixed bare coupling g with a value m such that 
X = holds up to a few sigma deviation. This is facilitated by the fact that rric 
seems to hardly vary with L at fixed g. For each of the observables considered, 

^The case u = v is special in the definition of S. It would depend on M{u) here but will not 
be required as we will not consider correlations with coinciding arguments (a; = 0). 
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for example those in ( HOI) , we also compute its derivative with respect to m. This 
is achieved by introducing another observable Wm, the derivative of the logarithm 
of the Boltzmann weight with respect to m, 



which is also easy to keep track of during the simulations. Now m-derivatives of 
not explicitly m-dependent observables O are given by connected correlations 

^{{0)) = {{OW^))-m){{Wr„)). (46) 

Thus we obtain estimates for dx/dm and, forming quotients, also dxn/dx (con- 
sidering g fixed). All these quantities are functions of primary expectation values 
and they are estimated including errors by the technique presented in [T2]. By 
exploiting the derivatives we can effectively 'post-run' vary m a little bit to solve 
the tuning problem of finding the value where x vanishes. This is equivalent 
to reweighting in m with a first order Taylor expansion of the Boltzmann factor. 
We thus for example obtain an estimate for = m + 6m with 

^dx/dm {{OWmSu,v)){{6u,v))-{{06u,.)){{Wm6u,v)) ^ ' 

with 

= n[2^--(o-o)] (48) 

a 

in this case. Since x is expected to be a smooth function of mL we expect that 
derivatives with respect to mL are order one. Then, for small x the next Taylor 
term can be roughly expected to be 0{x^^) which we monitor to be negligible 
compared to the statistical error of 5m and hence mc. Where this not true, the 
initial tuning of m was not good enough for first-order reweighting and the run 
has to be repeated closer to the chiral point. The range of possible corrections gets 
smaller as the number of lattice sites grows because of the variance of the estimate 
of the derivative. 

The entries in table [1] have been constructed in this way. Similarly as in the 
above discussion we have estimated the combination x„ — {dxn/dx)x to correct 
for tiny nonzero x- Often this systematic correction is insignificant, within one 
standard deviation, but the variance and thus the error is reduced due to anticor- 
related fluctuations of the terms. Each line in table [His based on a run of 8 x 10'' 
iterations, where an iteration performs elementary moves (counting u and v). 

The first series is at = 1, = 2 in the Thirring model. We expect no 
divergent coupling renormalization in this case and a continuum limit could be 
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m 


L 




Xi 


X2 


X3 


2 


1 


-0.4626 


8 


-0.46249(22) 


-0.0387(5) 


-0.0252(4) 


-0.0132(5) 


2 


1 


-0.4626 


16 


-0.46236(15) 


-0.0198(8) 


-0.0134(6) 


-0.0068(8) 


2 


1 


-0.4626 


32 


-0.46254(13) 


-0.0092(14) 


-0.0074(9) 


-0.0049(13) 


2 


1 


-0.4626 


64 


-0.46275(14) 

V / 


-0.0037(30) 


-0.0030(20) 


-0.0052(30) 


6 


0.3 


-0.5967 


8 


-0.59622(35) 


-0.0379(19) 


-0.0303(18) 


-0.0236(17) 


6 


0.3 


-0.5967 


16 


-0.59634(22) 


-0.0263(27) 


-0.0169(27) 


-0.0116(26) 


6 


0.3 


-0.5967 


32 


-0.59693(16) 


-0.0182(44) 

V / 


-0.0093(42) 


-0.0099(40) 


6 


0.3 


-0.5967 


64 


-0.59663(26) 


-0.015(25) 


-0.011(17) 


-0.012(11) 


6 


0.4 


-0.8007 


8 


-0.80075(36) 


-0.0558(19) 


-0.0389(18) 


-0.0286(17) 


6 


0.4 


-0.8007 


16 


-0.80060(21) 


-0.0312(29) 


-0.0200(28) 


-0.0121(26) 


6 


0.4 


-0.8007 


32 


-0.80077(15) 


-0.0246(47) 


-0.0119(45) 


-0.0014(42) 


6 


0.4 


-0.8007 


64 


-0.80081(17) 


-0.016(19) 


-0.015(14) 


-0.016(10) 


6 


0.5 


-1.0058 


8 


-1 00516f36) 


-0 0707fl9) 


-0 0423fl8) 


-0 0256fl8) 

KJ • KJ w \J KJ \ -L Ky t 


6 


0.5 


-1.0058 


16 


-1.00567(21) 


-0.0471(31) 


-0.0202(29) 


-0.0071(28) 


6 


0.5 


-1.0058 


32 


-1.00601(16) 


-0.0250(56) 


-0.0038(52) 


0.0066(50) 


6 


0.5 


-1.0063 


64 


-1.00612(33) 


0.059(69) 


0.050(37) 


0.025(20) 



Table 1: Simulation results for defined by x = and variant chiral order 
parameters Xn = 0(a) at m = m^. 



reached at fixed bare g, but the linearly divergent mass renormalization of Wilson 
fermions is present. As an asymptotically free example, data were taken for = 6. 
We notice in the table that the L dependence of rric is very weak. A systematic 
trend is barely visible. In figure [1] we plot our results for rric (taken from L = 64) 
for = 6, with the dominant one loop contribution subtracted, against g^. The 
dotted line corresponds to the curve 

3 

rric ~ ^m(")c/2", m^^^ = m^^^ m^^^ = -0.270, m^^^ = 0.189, (49) 

n=l 

where m^^ is the perturbative result fl66l) (see appendix A). 

It perfectly represents our data at the given error-level, but m^*) need not coincide 

with the true perturbative coefficients mc\i = 2, 3 that are not known to us. 

For N = 2 the only physical scale is L and we expect Xn to vanish proportional 
to a/L as chiral symmetry emerges in the continuum limit. A glance at the table 
shows that this is true within errors. For A^ ^ 3 asymptotic freedom holds and 
the situation is different. The continuum limit is expected to occur for g"^ \ 0. At 
least at large A^ [IJ we know that the physics is governed by an effective potential 
which for any g"^ has a critical Lc such that there is one minimum for L < L^. but 
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Figure 1: Values of the critical mass rric in lattice units at 
subtraction of the 1-loop contribution. 



6 after the 



two minima for larger L > Lc pL3j. At L = oo a fermion mass rriR is generated 
spontaneously and emerges as a nonperturbative scale, somewhat analogous to 
Aqcd- Therefore we expect for our A^ = 6 series a scaling behavior 

Xn = J^fnimRL) (50) 

with nontrivial but presumably smooth functions /„. 

At small g we expect ttirL to be tiny for L ^ 64 appearing in the table and 
we may replace the scaling function by /„(mijL) ^ /n(0) and thus the situation is 
similar to the Thirring model. This is still roughly so for the couplings ^ 0.5 
studied here. An investigation of the GN model at large A^ on the lattice with 
Wilson fermions has suggested that for N = 6, g"^ = 0.5 the value Lc may be not 
far from our L values, with quite some uncertainty however at finite A^. On the 
other hand the results in [13] suggest itirLc = 0(1) and we may be close to the 
onset of a change in the scaling behavior. 

In our data analysis we obtain information on autocorrelation times Tint.o for all 
primary and derived observables that we study. To avoid too much data handling 
we stored sub-histograms from 800 successive iterations resulting in a length of 10^ 
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of our time series. Between these blocked measurements autocorrelation times were 
small but sometimes still measurable. Under these conditions our error estimates 
are accurate at the percent level for the error of the error |l2j. To obtain more 
detailed dynamical information we made some shorter runs of 10^ blocks of only 
= 6 iterations each. This has yielded rint,x = 2.06(6), 3.47(14), 9.3(6), 69(10) 
for g"^ = 0.4 and L = 8, 16, 32, 64. The unit of Tint is now 'steps per site and 
flavor'. We see critical slowing down here, especially for the step 32 — > 64. This 
may be related to approaching a physically 'difficult' volume. The simulations of 
free fermions in |^7J have shown very little slowing down. The less favorable, but 
still not too bad, performance here is probably related to the flavors holding back 
each other. With growing coupling the different flavors get more more correlated, 
but they are only updated 'one at a time' with the others frozen until u = v is 
met. One could in principle try to influence the frequency of this by introducing 
a weight p{u — v) as done in [H] for the Ising model. In addition a second active 
flavor, needed anyway to simulate four-point functions, could be beneficial here. 

7 Conclusions and outlook 

We have successfully incorporated the GN interaction in the loop gas representa- 
tion of Majorana fermions and found the algorithm of |7j to remain rather efficient 
although some critical slowing down is visible. We have clarified the anomalous 
discrete chiral symmetry behavior in a finite volume. It can be used to define the 
chiral continuum limit with Wilson fermions by monitoring a topological observ- 
able X related to fiuctuating boundary conditions in the finite volume loop gas. It 
serves a purpose similar to the PCAC relation in the Schrodinger functional scheme 
for QCD and makes the required mass fine-tuning quite manageable. To complete 
the definition of a finite volume scheme for the GN model we at present still lack a 
good definition of a running coupling in a finite volume. The quantity proposed in 
[1] will presumably be a good choice but it requires the extension of the algorithm 
to also sample four-point functions. This seems feasible and will be required in any 
case to be able to probe mass ratios for which exact conjectures have been made 
|15] . It also seems advisable to implement on-shell Symanzik 0(a) improvement as 
in QCD. The action has to be augmented by only one term, for example, but 
in addition the elementary fermion operator ^ needs to be improved as it can mix 
with ^(^0- It would presumably be a good idea, interesting in itself, to first test 
this at large on the lattice with an analytic but nonperturbative calculation. In 
the end it may hopefully also be possible to shed some hght on the special cases 
= 3,4 which may be supersymmetric models in disguise |16j . 
From the algorithmic point of view we have seen here that with the sign problem 
absent in two dimensions (apart from a milder form in the finite size effects) we 
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could simulate a physically interesting interacting fermionic theory with the all- 
order strong coupling/hopping reformulation as an alternative to HMC methods 



A Some perturbative results 

In this appendix we explore in lattice perturbation theory to order g"^ some of the 
quantities that were studied numerically before. To evaluate quantities in leading 
order perturbation theory, it is useful to exploit the fact that the Grassmann 
integrals for each individual flavor are linear in each a{x) because it multiplies the 
bilinear at x whose square vanishes. Hence one derives 

-Z{a-z) = Z{a-z)K{a;z), ,:^^Z{a-z)=0 (51) 



da{x) da{xy 



with 



and 



k{a;z) = ^tTG{0;a;z) (52) 

G(x;a;.) = ^-"^^^^i"'^^^^^'"'^^ (53) 
Z{a;z) 

which is given in terms of propagators G{x; a; e) inverse to 'J^d^j, — \d^d*^ + cr with 
fixed boundary conditions e. Due to translation invariance we find 

^Z(m + ,a;.)|,.o = f^^(m;.) (54) 

with the volume V = LqLi, and similarly 

In addition we can make use of 

92 



;Z{a;z)G{x](T;z) = 0. (56) 



da{x) 

There is a small subtlety here. The vanishing of the doubly periodic partition 
function Z(0, (0, 0)) = trivially implies Z(0; z+) = Z(0; Z-). A singularity in the 
propagator G{x; m; (0, 0)) ~ (mLoLi)^^ + 0(m°) lifts however the zero-mode and 
yields an 0(a) difference 

2c 

lim [ir(m; z^) - k{m; z+)] = —=====. (57) 
—0 y/{Ll + Ll)/2 
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The constant derives from sums over momenta appropriate for the respective 
boundary conditions and the numerical Symanzik expansion (see appendix D of 
[n]) yields 

hm^^oZK (0,0))/(mL) 0.20668 , 

Coo = 7777: — ^ = 0.155609865 — h 0(L ). (58) 

L.^(o,o)^(0;^) L 

Here and below we perform the numerical evaluations for Lq = Li = L. If we 
write K{0; z_) — K{0; z+) below we refer to fl37|) . 

A.l Chiral order parameter x 

With the help of the above formulas it is now rather easy to derive 

Z^^(m, g) = [Z{m- z^)f |l + ^iV(iV - l)Vk{m- z^f^ + 0{g^) (59) 
and hence 

X = ln[Z(m; z.)/Z{m- z+)] + ^(iV - l)V[k{m- z^f - k{m- z+f] + 0(/) (60) 

The perturbative Ansatz for the critical mass of Wilson fermions reads 

mc(^^) = c/^mW + g^mf^ + .... (61) 
Here we may use to leading order 

Zf'{g^m^\Q) = [Z(0;z±)]^{l + /m«Wir(0;^±)}. (62) 

With X vanishing at tree level we now derive (up to 0{g^)) 

X = ^ [ir(0; z^) - kiO; z^)] {{N - 1) [i^(0; z^) + 7^(0; z+)] + 2m« } . (63) 

The condition that x vanishes as a chiral order parameter fixes the one loop mass 
renormalization to 

m^c^ = -^^[k{0; z.) + k{0; z+)]. (64) 
Using the definitions this may also be written as 

,r,_ iv-i S.^w)g(0;^)fG(0;0;^) 

' " ~ E.,<„,„,Z(0;.) '''' 

and the numerical value is 

2439 

m(^) = -(N - 1)0.3849001795 - — — + 0(L-^). (66) 
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A. 2 Two- point correlator 

For the correlator ( l24l) (for one flavor a = b) the following simple result is obtained 
(e(x)e(O)). = ^ |l + 9\N - l)k{m- z) — I G{x- m- z) + Q{g'). (67) 

In this formula the partition function ratio can also be expanded and contributes 
in general another 0{g^) term. At the chiral point however it is unity for z = 
and trivially so for z+. The correlations defined in this paper now read 

kvixo) = |l + g'{N - l)k{m; ^+)^} kP{xo) + Oig') (68) 

and 

ks{xo) = \ J2 e^^^--|l+/(iV-l)ir(m;^,)A|40,.)^^^)_^Q(^4)_ ^gg) 

T = +,- ^ ^ 

Using results from appendix A in [7j we find (for < Xq < Lq) 



= -xi Z(m,(l,0)) 



cosh[c<j(Lo/2 — xq)] 
1 + m)Z{m; z+) " P "'^ cosh[cuLo/2] 

-Z(m,(0,0)) ""^["[tf 7r°H (70) 
^ ' sinh[cuLo/2] / ^ ' 

where u = ln(l + m) is the pole mass. For small m we obtain 

(1 + m)k^y\xo) ~ 2cio(l - coomLo) - 2coom{Lo - 2xo) + O(m^) (71) 

with another constant 

^w = ^ —— = 0.313557560 — h 0(L ). (72) 

E.^(o,o)^(0;^) 

For the scalar correlator we derive 

kf^\xo) = 1 x|z(m,(l,0))^^"^["^^°/'-^°)] 



{l + m)Z{m;z±) { ' ' cosh[c^Lo/2] 

T^(m,(0,0)) ^"^^["ftf (73) 
^ 'V ' sinh[tuLo/2] J ^ ^ 

In this case we find for small m 

(1 + m)kf^\xo) = 2ciom(Lo/2 - a;o) T 4coo + 0{m^). (74) 
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If we combine the above results we find that the alternative order parameters 
Xn vanish also at 0{g'^) if we set m = g'^rni \ i.e. the 0(a) contributions seen 
numerically start (presumably) at O ((?"') only. 
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